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Abstract 

We simulate 50 off-lattice DLA clusters, one million particles each. Lhe 
probability distribution of the angle of attachment of arriving particles with 
respect to the local radial direction is obtained numerically. For increasing 
cluster size, N, the distribution crosses over extremely accurately to a co- 
sine, whose amplitude decreases towards zero as a power-law in N. From 
this viewpoint, asymptotically large DLA clusters are locally isotropic. This 
contradicts previous conclusions drawn from density-density correlation mea- 
surements [P. Meakin, and T. Viscek, Phys. Rev. A 32, 685 (1985)]. We 
present an intuitive phenomenological model random process for our numeri- 
cal findings. 
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I. INTRODUCTION 



The diffusion limited aggregation (DLA) M is a fractal growth model exhibiting great 
complexity 0. The properties of asymptotically large clusters are of fundamental impor- 
tance. Some of these properties are independent of many details of the aggregation rules 
and are expected to be shared by large clusters grown in various experimental or natural 
situations. Efficient algorithms and improved computational facilities has enabled the gen- 
eration of very large off-lattice clusters with more than 100 million particles ||. However, 
many asymptotic properties of DLA are still unclear. There have been many debates on 
basic issues including the scaling behavior of large clusters and the multi-fractal properties 
of the growth probability measure 0. 

This work concentrates on one of the most fundamental geometrical aspects, namely the 
local isotropy, of large off-lattice DLA. Meakin and Viscek [|J found that the two points 
density-density correlation is anisotropic. The tangential correlation decays algebraically 
with exponent a± ~ 0.41, which is different from the radial exponent an ~ 0.29. As a re- 
sult, the geometry, in particular the orientation, of a segment of the cluster is related to its 
position inside the cluster. This behavior does not hold for some deterministic fractal curves 
containing spirals such as the Koch curve || in the limit of small segments. Asymptotically 
large DLA would be locally isotropic only if the two very different exponents converge to the 
same value. Such a trend was not identified f| and the anisotropy has been assumed to per- 
sist asymptotically flPH. Here, we apply a substantially more accurate method and obtain 
the first unambiguously evidence of a systematic decrease in the local anisotropy of DLA 
for increasing cluster size. The trend of the decrease strongly indicates that asymptotically 
large clusters are locally isotropic, contradictory to previous conclusions (|J. 

Specifically, we investigate the orientation of particle attachment. For each new particle 
in the cluster, let R' and R be the position vectors relative to the center the cluster of 
respectively the new particle and its parent, where the parent is the particle in the aggregate 
upon which the new attachment is made. We define the center of the cluster to be at the seed. 
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The vector R characterizes the local radial direction of this sticking event and r = R! — R 
gives the sticking direction. The angle of attachment 9 (—71 < 9 < tt) is defined as the 
angle measured counter-clockwise from R to r. Attachment with with | 9 \ < tc/2 can 
be described as forward. We focus on the probability distribution P(9, N) of 9, which, in 
general, depends on the number of particles, N, in the cluster. If DLA were compact with 
nearly circular boundary, the above definition implies that 9 is always close to zero and a 
backward attachment is geometrically impossible. The distribution P(9, N) should peak 
at 9 = 0. Although, DLA is far from being compact, at least for small clusters, forward 
attachment is favored for similar geometrical reason. 

II. SIMULATIONS AND RESULTS 

We compute the distribution P(9, N) for N = 150 x 2 k , where k=0 to 12 corresponding 
to N=150 to 614400. For each k, we histogram the values of the sticking angle 9 for the 
(100 x 2 k + l)-th to the (200 x 2 fc )-th particles. Proper normalization gives P(9,N). The 
symmetry P(9, N) = P(—9, N) is used. We averaged the results over 50 clusters. The 
statistical error is estimated from the sample to sample fluctuations. 

Figure [l] shows P(9, N) for 5 values of N. We used 16 bins to histogram. The error bars 
are smaller than the symbols except for N = 150. The errors are as small as about 0.2% at 
N = 614400, where our data is most accurate. For small N, forward sticking dominates so 
that P(9, N) peaks at 9 = as expected. Backward stickings are rare. The peak broadens 
as N increases. At N = 614400, the ratio P(±7r, N)/P(0, N) between perfectly backward 
and forward attachments becomes about 0.31. Figure [I] show least square fits to: 

P(9, N) = l/2?r + ai(N) cos(9) (1) 

for the 3 larger values of N, where 1/27T ensures normalization. At N = 614400, the quality 
of the fit amazing. There is apparently no systematic trend of deviation, as the data points 
scatter around the fitted curve by amounts comparable to the statistical errors. The fit is 
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still good at N = 76800, but at iV = 9600, there are noticeable systematic deviations. The 
quality of the fit is restored when we consider one more term in the cosine expansion: 

1 oo 

P{6, N) = — + J2 a n (N) cos{ne) (2) 

Z?l n=l 

Figure |I] shows the two parameters fits for N = 150 and 1200. More generally, we fitted our 
complete data set of P(9, N) with the first 4 cosine terms. Figure 0(a) plots the values of 
the amplitudes a n (N) in semi- log. We used histograms of 64 bins and the result agrees with 
those with 16 bins. We also computed a n (N) by Fourier transform and obtained again the 
same results. Figure 0(b) and (c) shows respectively a%(N) and a,2(N) in log- log plots. The 
linearity at iV > 300 implies: 

a n (N)^A n N~^ (3) 

for n = 1 and 2, where 71 = 0.0997(3), A x = 0.309(2), 72 = 0.67(3) and A 2 = 1.2(3). The 
bracketed values are the fitting errors. For n > 3, the measured a n (N) is not precise enough 
for a test of the above algebraic decay. 

We also studied the closely related problem of the branch orientation of DLA. We adopt 
the Horton-Strahler's scheme of branch ordering |7J. The smallest branches without side- 
branch are assigned order 1 and the main stems have the highest orders. Other authors have 
used slightly different schemes |||| , which should give similar results to those reported here. 
From visual examinations of the branches shown in Ref. ||, it is evident that those with 
the highest order branches are nearly radial and the directedness decreases for lower orders. 
Quantitatively, we define the branch orientation angle <fi to be measured counter-clockwise 
from the position vector of the base of the branch to the branch orientation vector pointing 
from the base to the tip. Figure ^| plots the probability distribution P(0, s) of </> for the 
branch order s — 1 to 5. The data was averaged over 50 clusters, one million particles each. 
The distribution P(<p, s) can be approximated reasonably well by a cosine curve for n close 
to 1. In the cosine expansion of P(<f>, s), the coefficient of the cos((p) term decreases rather 
quickly with s, while that of cos (2<j)) is negative and rising up towards 0. The coefficients 
as functions of s are not well described by simple functional forms. 
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III. PHENOMENOLOGICAL MODEL 



We now propose an intuitive phenomenological picture for the properties of the angle of 
attachment. For convenience in presentation, we examine the following approach of DLA 
growth, instead of the more efficient random walker method. To add a particle to the 
cluster, the Laplacian potential is solved to obtain the particle probability flux lines. The 
new particle is launched randomly with uniform probability on a big circle inscribing and 
far away from the cluster. It subsequently traces the flux line deterministically towards the 
cluster until it hits and become part of it. In this scheme, the direction of the attachment 
is simply the tangent to the flux line at the point of contact. Figure § shows 200 flux lines 
with uniformly spaced starting positions far away. They are all equally likely to represent 
the next growth step. We will concentrate on the geometry of these typical flux lines which 
are relevant to the growth. Lines with very low probability weight, which can have very 
different geometry |l(J, are neglected. 

The flux lines proceed nearly exactly radially inwards until they are close enough to the 
cluster to be influenced by the geometry of the main branches. At this stage, the details of 
the side-branches are still unimportant. Depending on the position of the flux line relative 
to the nearest main branch, it can keep proceeding radially towards the tip or make a turn 
to approach the branch from either side. After advancing further, the geometry of the 
nearest side-branches becomes important. Again, the flux line approaches either towards 
the tip or one of the sides of the side-branch. Similar situation recurs for side-branches of 
the side-branches until the geometry of individual particles becomes relevant. 

The typical total number of turns thus equals the number of levels of side-branching, 
which is proportional to r = In N The angle of attachment is approximately the sum 

of the angles of turning. Suppose that when r is increased by At, there are, on average, one 
more turn in the flux lines. We approximate the corresponding evolution of the distribution 
of the angle of attachment by: 




(4) 
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where there is a reparametrization with r. Implicit assumptions includes the uncorrelation of 
the turns and their isotropy so that the propagator G only depends on 9 — 9'. Unfortunately, 
a direct numerical verification of Eq. (|4]) with the determination of G requires computation 
of the flux lines for very large clusters, which is well beyond the capability of the relax- 
ation method 0. The validity of the assumptions can only be justified by comparing their 
predictions with simulations. 

We consider first the simplest case where G is a Gaussian centered at A9 = with width 
smaller than tt. Equation (f|) reduces to: 

§;P{0,T)=u£ s P{0,T) (5) 

where v is the diffusion constant. The set of Eqs (0) and (|3]) is indeed a solution to Eq. 
(H). They predict 7„ = vn 2 and thus the amplitudes of the higher harmonics decay much 
faster. This is consistent with the practically zero values of the measured amplitudes for 
n > 3. Moreover, the result 72/71 = 4 is in the same order as the numerical value 6.7(3). 
For general G, we expand: 

1 00 G 

G{A9) = — + y ^cos{nA9) (6) 

27F n=l 71 

Substitution into Eqs (§) and (§) gives 7„ = — lnG n /Ar. Consider the next simplest form 
besides the Gaussian: G(A9) ~ exp[— (A9/9 ) 4 ]. By matching 71 and 7 2 with the numerical 
values, we obtain respectively 9q = 94° and At = In 130. It means that a flux line make a 
turn of typically not more than 94° per increase in iV by a factor of 130. For a cluster of 
N = 10 6 , each flux lines has roughly N t ~ 2.8 turns. We also tried G in a very different 
form of the sum of two Gaussian peaks. It gives Nt — 4. 

Recall that the number of turns N t equals approximately the number of levels of side- 
branching probed by the flux lines. The maximum branch order in a cluster of a million 
particles is typically 9 [[| . However, most flux lines terminate around the tips of the cluster 
so that the higher order branches, which reside close to the center, are irrelevant. As a result, 
Nt 9. It seems that N t ~ 3 or 4 obtained from the very crude models are reasonable 
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approximates. Even for large clusters, N t is very small. This is the cause of the extremely 
slow evolution of P(9, N). 

The length scale of the problem decreases as we consider deeper levels of side-branching. 
In our construction with the modified Gaussian, the displacement of the flux line in between 
turns typically decreases by a factor of 130 1 / 13 ~ 17 after every turn, where D ~ 1.72 is the 
fractal dimension of DLA ||. every flux line converges rapidly to its point of contact. For 
extremely large clusters and length scales in between that of the cluster and the individual 
particles, segments of flux lines containing their points of contact are statistically similar 
to each other after proper rescaling. This scaling property is a consequence of the approxi- 
mation that the propagator G is independent of N, based on the assumption that the lines 
progress in self-similar environments. Dependence of G on N should reveal in discrepancies 
from Eq. (0), which corresponds to any non-linearity in the log-log plot in Fig. |2](a) and 
(b). Within our accuracy, we cannot identify any unambiguous non-linearity for iV > 300. 
Our measurement is not sensitive to the well known deviations from self-similarity of DLA 

00- 

Although the uncorrelation of the turns seems to be a reasonable approximation, the 
are obviously strong correlations within a turn so that Eq. (^) cannot be reformulated with 
smaller At corresponding to fractions of a turn. Indeed, G(A8, At/2) is not well defined in 
both models with the modified Gaussian and the double Gaussian, since G3 < 0. 

The behavior of the distribution P(<f), s) of the branch orientations can be accounted for 
similarly. The main branches with the highest orders are nearly radial. Their side-branches 
are attached to either sides with a subtended angle around ±40° |§. This contributes to 
a step in the randomization of the orientation analogous to the turns of the flux lines. 
Therefore, much of above discussion applies. However, the correlation in the randomization 
process is stronger here, since a branch physically excludes some of the side-branches of 
smaller orders to have the same orientation. It is the exclusion by the radial main branches 
which causes the dips of P(4>, s) at = in Fig. [| This correlation also lead to the break 
down of any analogous power law decay of the amplitudes in the cosine expansion of P(<f), n). 
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IV. DISCUSSIONS 



The local anisotropy of DLA is far from being stabilized for clusters of size N < 10 6 . 
Assuming that Eq. (Q) can be extrapolated to iV — > oo, we get lim^v^oo P(Q, N) = l/2n and 
DLA is locally isotropic asymptotically The reliability of this extrapolation deserves special 
attention due to its simplicity, excellent agreement with simulation for 300 < N < 10 6 , 
and consistency with an intuitive phenomenological model. At finite N, the anisotropy 
can quantitatively be expressed by the proportion, Pp, of forward sticking events given by 
P F ~ 0.5 + 2a 1 N~' fl from Eq. (g). For N = 10 3 to 10 6 , P F decreases from 80% to 65%. 
The anisotropy is still strong and agrees with the density correlation measurements [[§]. 
Achieving approximate isotropy with Pp ^5 55% requires N > 10 11 ! 

Our results on P(6, N) is one of the most accurate non-trivial measurements ever done on 
DLA. In particular, the amazing accuracy in the one parameter fit to the simple analytic form 
in Eq. ([!]) at iV = 614400 is rare. The cos(8) term, which is the dominant term representing 
the anisotropy, was not identified before because forward and backward attachments are not 
properly discriminated in the correlation measurements 

We hope to motivate similar measurements for diffusion limited deposition (DLP). If 
analogous anisotropy tends to zero asymptotically, very large DLP might not be self-affine 
\T2\. The same approach might also be attempted for the dielectric breakdown model 0. 

The numerical part of this work was done with Henry Kaufman and Benoit B. Mandel- 
brot, who are gratefully acknowledged. 
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FIGURES 

FIG. 1. Probability distribution P(9, N) of the angle of attachment for different cluster sizes 
N. Also shown are one parameter cosine fits for the data at N = 9600, 76800 and 614400 and two 
parameters fits for N = 150 and 1200. 

FIG. 2. Amplitudes a n (N) of the cosine terms in the expansion of P(9,N) as a function of 
cluster size N: (a) a\ to 04 in semi-log plot; (b) a\ and (c) 02 in log-log plots respectively. 

FIG. 3. Probability distribution P(4>, s) of the branch orientation angle 4> for branches of order 
s and fits with four cosine terms. 

FIG. 4. Example of 200 particle probability flux lines with uniformly spaced starting position 
on a circle of radius ARq, where Rg is radius of gyration of the one million particles DLA cluster. 
The potential was solved by over-relaxation method in a circular region of radius ARq on a grid 
with lattice spacing 8Rg/7000. Displayed region has radius 2Rq- 
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